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THE  DYNAMIC  BEHAVIOR  OF  WATER  VAPOR 
ON  ACTIVATED  CARBON 


1.  INTRODUCTION 

Water  adsorption  on  activated  carbon  is  of  particular  interest  for  air 
purification  applications  where  contaminants  must  be  removed  from  humid  air. 
Water  vapor  adsorption  can  limit  filter  performance  against  challenges  of  light 
gases  by  displacing  previously  adsorbed  contaminants  or  by  reducing  the  available 
adsorption  space.  A  complete  understanding  of  the  single  component  adsorption 
behavior  of  water  vapor  on  activated  carbon  can  serve  as  a  basis  to  develop  a 
theory  for  multicomponent  behavior^  This  study  considers  in  detail  the  adsorption 
of  water-on  BPITactivated  carbon. 

The  adsorption  behavior  of  any  vapor  on  an  adsorbent  can  be  understood 
through  the  adsorption  equilibrium  relationship  and  the  material  and  energy 
balances.  Water  adsorption  on  activated  carbon  exhibits  a  pronounced  hysteresis 
that  complicates  the  analysis  of  desorption.  Doong  and  Yang1  present  a  model  for 
multicomponent  adsorption  of  water  and  hydrocarbons  on  activated  carbon 
including  equilibrium  and  adsorber  dynamics.  Their  work,  however,  neglects  the 
energy  balance  and  does  not  present  data  on  breakthrough  runs.  The  purpose  of 
this  paper  is  to  develop  and  solve  the  material  and  energy  balance  equations  so  the 
results  can  be  compared  to  actual  breaktlirough  experiments.  The  results  presented 
here  for  adsorption  dynamics  can  be  extended  to  desorption  dynamics.  Two 
approaches  to  analyzing  the  water  adsorption  process  on  activated  carbons  are 
considered.  These  are  the  stage  model  approach  and  the  equilibrium  theory 
approach.  In  both  cases  a  mathematical  model  is  used  to  describe  the  dynamic 
behavior  of  water  vapor  challenging  a  bed  of  activated  carbon.  Three  breakthrough 
experiments  are  performed  that  result  in  three  different  types  of  waves.  The 
characteristics  of  each  wave  type  and  their  impact  on  filter  performance  are 
discussed. 

2.  THEORY 

2.1  Material  and  Energy  Balances. 

The  equations  to  describe  water  adsorption  on  activated  carbon  should 
include  material  and  energy  balances  as  well  as  the  equilibrium  relationship.  The 
easiest  way  to  develop  these  equations  is  to  describe  the  bed  as  consisting  of  a  solid 
phase  and  a  vapor  phase.  The  stationary  (solid)  phase  consists  of  the  solid  porous 
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particles  and  the  adsorbate.  The  vapor  phase  consists  of  the  inert  carrier  gas  and 
the  adsorbable  vapor  in  both  the  intersticies  of  the  bed  and  the  pores  of  the 
particles. 


The  governing  differential  equations  for  the  conservation  of  mass  and 
energy  respectively,  neglecting  dispersion  and  channeling,  of  a  fixed-bed  adsorber 
can  be  written  as  follows: 


dvc 

dz 


+  € 


dc 

dt 


+  Pb 


ia.  _ 


at 


=  o 


(i) 


d(vPfhf)  d(pfhf) 


dz 


+  fi¬ 


at 


+  Pb' 


djh 

at 


=  -U  AT 


(2) 


The  fluid  phase  and  stationary  phase  enthalpies  are  defined  by  : 


hf  —  cf(T  -  Tref) 


(3) 


—  (cs  +  ca<l)(T  -  Tref)  —  /o  A  dq 


(«) 


Here  c  is  the  vapor  concentration,  q,  the  adsorbed  phase  concentration  in  moles  of 
adsorbate  per  unit  mass  of  adsorbent,  T,  the  temperature,  z,  the  distance  down  the 
tube  length,  t,  the  time,  e,  the  void  fraction  of  the  bed,  and  U,  the  overall  heat 
transfer  coefficient  of  the  filter  tube. 

The  rate  at  which  material  is  transferred  from  the  passing  vapor  stream  to 
the  surface  of  the  carbon  is  governed  by  a  film  mass  transfer  resistance  which  can 
be  written  as  follows: 


Pb  =  (M)  (c  -  (5) 

Equation  5  assumes  Glueckaufs  linear  driving  force,  where  c*  is  the  vapor  phase 
concentration  at  the  surface  of  the  particle. 

Dimensionless  variables  can  be  introduced  to  simplify  equations 

1,  2,  5: 
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T  = 


tVpt 

L 


(6,7,8) 


U' 


2  U  L 
c  v0  rj 


(9) 


(kva)' 


(10) 


The  choice  of  the  dimensionless  parameter  r  has  special  physical 
significance.  It  represents  the  number  of  superficial  column  volumes  that  have  been 
fed  to  the  bed  at  any  given  time. 

In  evaluating  solutions,  parameters  appearing  in  equations  1,  2,  and  5  can 
be  treated  as  functions  of  concentrr  Lion  and  temperature.  However,  to  simplify  the 
development,  several  assumptions  can  be  made,  namely: 

(a)  ail  heat  capacities  are  constant  with  mean  values  taken. 

(b)  the  vapor  phase  is  ideal  and  the  effect  of  composition  on  the  vapor- 
phase  density  and  velocity  is  small. 

(c)  the  adsorption  process  considered  here  occurs  at  low  gas  phase 
concentrations  so  that  the  effect  on  velocity  and  density  is  negligible. 

Thus,  the  dimensionless  velocity  and  the  fluid  density  can  be  written  as: 


* 

V  = 


(11,12) 


For  adsorption  processes  involving  gases  at  low  to  moderate  pressures 
(which  would  be  typical  of  all  air  purification  applications)  the  rate  of  accumulation 
in  the  vapor  phase  is  small  compared  to  the  accumulation  in  the  solid  phase.  This 
leads  to  the  elimination  of  the  second  term  in  equations  1  and  2.  The  system  of 
governing  differential  equations  may  then  be  expressed  as: 


d(v*c) 


+  Pb 


dr 


=  0 


(13) 
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d(vVffy)  dhg 

~~aT~  ^  IT 


-U'AT 


(14) 


Pb^  =  (kva)'  (c  —  c*>  (15) 

Equations  13,  14,  and  15  are  three  independent  equations  with  four 
dependent  variables  (c,  c  ,  q,  T).  Another  equation  is  required  in  order  to  specify 
the  system,  namely  the  adsorption  equilibrium  relationship: 


q  =  f(c*,T) 


(16) 


The  solution  to  equations  13  -  15  with  16  requires  the  imposition  of 
appropriate  initial  and  boundary  conditions.  For  systems  considered  here  the.  initial 
bed  state  and  the  feed  concentration  are  constant  throughout  the  experiment. 
Mathematically  this  can  be  expressed  as: 

c  =  c0  and  T  —  T0  at  r  =  0  for  0  <  f  <  1 
c  —  cfeed  T  =  Tf^  at  f  =  0  for  all  r  >  0 

Two  approaches  are  taken  to  solve  the  governing  set  of  differential 
equations,  13  - 15,  with  16:  the  stage  model  and  equilibrium  theory.  The  stage 
model  is  more  flexible  and  can  easily  treat  systems  with  non-constant  initial  and 
boundary  conditions.  Generally,  equilibrium  theory  is  more  insightful  about  the 
nature  of  the  adsorption  waves,  however,  it  cannot  be  used  where  mass  and  heat 
transfer  effects  must  be  considered.  Both  techniques  are  presented  here  in  order  to 
gain  a  complete  understanding  of  the  process. 

2.2  Stage  Model  Approach. 

Equations  13  -  15  when  coupled  with  the  equilibrium  relationship  16  can 
be  reduced  to  a  set  of  ordinary  differential  equations  using  a  stage  model 
approach.’3  Each  stage  is  considered  to  be  well  mixed  with  the  effluent 
concentration  and  temperature  of  stage  i  used  as  the  feed  concentration  of  stage 
i+1.  Friday  and  LeVan  have  reported  a  successful  method  for  reducing  the  partial 
differential  equations  13  -  15  to  ordinary  differential  equations  by  writing  the  spatial 
derivatives  in  equations  13  and  14  in  finite  difference  form  using  backward 
differences. 
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dqi 

dr 


1 


(18) 


dTj 

dr 


1 


Pfocf 


Pb(cS+caqi)l  Af 


(Ti-1  -  TO 


do;  i 

-PbMTi-T^J-Aij-^-U'AT 


(19) 


The  resulting  equations  18  and  19  axe  solved  using  a  Runge-Kutta  4th  order 
integration  algorithm  with  twenty  stages. 

2.3  Equilibrium  Theory. 

The  development  of  the  equilibrium  theory  solution  requires  the  mass  and 
heat  transfer  resistances  to  be  negligible.  Therefore,  the  driving  force  terms  on  the 
right  side  of  equations  14  and  15  are  set  to  zero. 

The  equilibrium  theory  solution  to  13,  14  and  16  requires  the  use  of  the 
method  of  characteristics  that  has  been  presented  by  numerous  authors.5-10  If 
concentration  of  the  material  in  the  bed  can  be  expressed  as  a  function  of  distance 
and  time  by, 


c  =  fU,r)  (20) 

then  the  total  differential  of  c  can  be  written  as: 

<21> 

The  solution  along  regions  of  constant  concentration  require  that  dc=0  so  that 
equation  19  can  be  rearranged. 


_d£  _  dc/dr 
dr  dc/df 


(22) 


Now  the  material  balance  equation  13  can  be  rewritten  when  the  equilibrium 
relationship  16  is  substituted  for  the  partial  of  q  with  respect  to  time. 
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(23> 


—  + 


Ph 


Ai. 


d(v*c) 


dc 

3t 


\ 

Equation  23  can  be  rearranged  into  a  form  compatible  with  the  right  hand  side  of 
equation  22. 


dc/dr  _  _1_  d(v*c) 
dc/df  pb  dq 


(24) 


When  equation  24  is  substituted  into  equation  22  the  following  expression  for  the 
characteristic  direction  is  obtained. 

.L  d(v*c) 
dr  Pb  dq 


The  left  hand  side  of  equation  25  represents  the  speed  with  which  a  given 
concentration  travels  down  the  bed,  while  the  right  hand  side  of  equation  25  is 
the  slope  of  the  equilibrium  relationship. 

The  energy  balance  when  written  without  the  heat  transfer  resistance 
takes  the  following  form: 


d(vVfhf) 

3?  ' 


+  Pb' 


d_K 

8t 


(26) 


Using  the  same  procedure  as  for  the  material  balance,  the  characteristic 
equation  takes  the  following  form: 


_d£  _  _1_  d(v*pfhf) 
dr  pt  dhg 


The  material  and  energy  balances  can  be  combined  for  the  initial  value 
problem  to  provide: 
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L  _  dq  _  dhg 
C  ^  d(v*c)  ^  d(v*p(h{) 


(28)  ' 

\ 

! 


The  two  enthalpy  terms  on  the  right  of  equation  28  can  be  written  as 
functions  of  q  and  T.  Thus  equation  28  can  be  rearranged  in  terms  of  the 
directional  derivative  dq/dT  to  take  the  following  form: 


D(v*c)  dhg 

’  > 

i - 

a. 

.o 

1 _ 

dT 

«.  y 

d(\  c)  dkg  d(v*c)  dbg 
dq  ST  +  ST  Sq 


d(v*Pfhf) 

dT 


dq  S(v*c)  dh3 
dT  +  ST  ST 


0 


(29) 


Equation  29  is  a  quadratic  in  dq/dT.  The  two  solutions  correspond  to  the 
characteristics  emanating  from  a  given  point  of  c,  q,  and  T.  By  starting  from  the 
initial  and  feed  conditions  the  solution  to  equation  29  is  obtained.  The  resulting 
slopes  are  then  plotted  in  the  hodograph  planes  of  q  versus  T  and  the  q  versus  v  c 
until  the  characteristics  intersect.  Implementation  of  equilibrium  theory  to  obtain 
the  breakthrough  curves  requires  adherence  to  two  rules.  A  wave  must  be  gradual 
if  possible.  Slopes  in  the  v  c  plane  must  increase  from  initial  bed  condition  to  feed 
condition.  If  the  slope  of  the  characteristic  does  not  increase  continuously  in  going 
from  initial  to  feed  conditions  then  the  larger  concentrations  would  overrun  the 
smaller  concentrations.  This  is  an  impossible  situation.  Earlier  studies  by 
Devault11  and  Lax12  showed  that  these  conditions  result  in  the  formation  of  a 
shock.  The  shock  path  must  satisfy  the  material  and  energy  balances  as  given  by 
the  following  relationship: 


L  -  Aq  _  Ahg 

£  ^  A(v*c)  ^  A(v*pfhf) 


(30) 


This  is  similar  to  equation  28  except  that  here  a  chord,  rather  than  the  slope  of  the 
characteristic,  is  used  to  determine  the  speed  of  the  wave  in  the  bed. 


3.  MATERIALS  AND  METHODS 


The  apparatus  used  for  the  breakthrough  experiments  is  shown  in 
Figure  1.  Air,  used  as  the  carrier  gas,  passes  through  a  molecular  sieve  dryer  and  a 
carbon  Filter  to  remove  water  and  organics.  The  main  metering  valve  regulates  the 
flow.  After  the  clean  air  passes  through  the  flowmeter,  a  portion  of  the  stream  is 
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saturated  in  a  water  bath.  The  desired  relative  humidity  is  achieved  in  the 
following  manner.  The  metering  valve  and  the  bypass  valve  are  set  such  that  the 
ratio  of  the  their  flows  produces  a  relative  humidity  near  the  desired  value.  The 
humidity  analyzer  and  controller  (EG&G  911  Digital  Humidity  Analyzer)  measures 
the  relative  humidity  of  the  feed  stream  based  on  the  temperature  of  the  feed 
stream.  If  the  relative  humidity  is  below  the  setpoint,  the  water  bath  heater  is 
activated.  Since,  it  is  a  simple  on-off  controller  with  no  cooling  action,  the  initial 
setting  of  the  bypass  flow  must  produce  a  lower  than  desired  RH.  The  humidified 
air  flows  through  the  carbon  bed  and  the  effluent  relative  humidity  is  measured  by 
a  second  EG&G  humidity  analyzer.  Periodically  during  each  breakthrough  run  the 
bed  was  removed  from  the  flowing  air  stream  in  order  to  weigh  it. 

The  carbon  sample  used  was  12x30  mesh  BPL  Lot  #7502.  Calgon 
literature  states  that  the  bulk  density  of  the  carbon  is  480  kg/m3.  The  particle 
diameter  is  taken  to  be  0.001  m.  The  sample  for  each  experiment  is  dried  for  at 
least  two  days  at  110  °C.  The  carbon  is  then  poured  into  a  0.0415  m  diameter 
plexiglass  tube  to  give  a  bed  depth  of  approximately  4  cm.  The  overall  heat 
transfer  coefficient  of  the  plexiglass  is  taken  to  be  0.093  k J /m2  K  s. 

4.  RESULTS  AND  DISCUSSION 

4.1  Water  On  Activated  Carbon. 

The  breakthrough  concentration  and  temperature  have  been  measured 
using  the  following  feed  and  initial  bed  conditions:  50%  RH  feed  to  a  clean  bed,  80% 
RH  feed  to  a  bed  equilibrated  at  50%  RH,  80%  RH  feed  to  a  clean  bed.  These 
conditions  have  been  chosen  to  illustrate  the  three  different  wave  types  that  can  be 
formed  in  an  adsorption  bed.  Rhee  et  al.  describe  the  three  wave  types  as- 
gradual,  abrupt  and  combined  in  their  equilibrium  theory  analysis.  The  important 
features  of  each  are  discussed  as  each  experiment  is  analyzed.  Table  1  is  a  summary 
of  the  operating  conditions  for  each  experiment.  Table  2  lists  the  heat  capacity 
values  used  in  this  study. 
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TABLE  1.  Breakthrough  Experiment  Conditions 


Experiment 

Feed 

Initial  Bed  Conditions 

V 

T 

amb 

RH% 

T(K) 

RH% 

T(K) 

Carbon  (g) 

(m/sec) 

(K) 

1 

50 

298 

0 

297.5 

23.96 

.0616 

297.5 

2 

80 

297 

50 

297 

23.96 

.0616 

297 

3 

80 

298 

0 

298 

25.85 

.123 

298 

TABLE  2.  Heat  Capacities 


c 

a 


c 

8 


=  1.04  kJ/kg  K 
=  2.7  kJ/kg  K 
=  0.075  kJ/kg  K 
=  298°  K 


In  order  to  implement  the  stage  model  and  the  equilibrium  theory  model 
the  form  of  the  equilibrium  relationship  must  be  established.  It  is  proposed  that  a 
modification  to  the  Hacskaylo-LeVan  isotherm13  would  give  a  reasonable  fit  to  the 
data.  The  form  of  the  equation  is  as  follows: 


In  (p) 


A'  - 


B' 

C'  +  T 


(31) 


where 

A'  =  A  +  (K10  +  l)ln(0)  B'  =  B  C'  =  C  +  K2(l  -  6)  (32) 

and 
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(33) 


5  =  -i- 

9sat 


A  -  RT2 


d[n  (eI 

dT 


=  RT2- 


B' 


(C'  +  T)2 


(34) 


The  following  values  axe  used  for  water  adsorption  on  BPL: 

A  =  18.3036  B  =  3816.44  C  =  -46.13 

K1  =  2.56  K2  =  37.9  q^  =  0.0217 

The  Antoine  parameters,  taken  from  Reid  et  al.,14  give  the  vapor  pressure 
in  millimeters  of  mercury.  The  isotherms  of  water  on  BPL  have  been  measured 
previously  by  Mahle  and  Friday15  over  the  temperature  range  15  -  35  °G.  Only 
adsorption  data  is  considered  here;  desorption  is  not  considered.  The  result  of 
correlating  the  adsorption  data  Vising  equations  31  -  33  is  presented  in  Figures  2  and 
3  in  the  modified  isotherm  plane  of  q  versus  v  c.  Figure  3  is  an  enlargement  of 
Figure  2  in  the  low  loading  region. 

The  value  of  ky,  the  mass  transfer  coefficient,  can  be  calculated  using  the 
Sherwood  number  correlation  of  Ranz  and  Marshall.16  Calgon  literature  indicates 
that  an  appropriate  value  for  a,  the  surface  to  volume  ratio,  is  3400  m'1. 

4.2  Breakthrough  of  50%  RH  Feed  To  A  Clean  Bed. 

The  breakthrough  results  for  a  challenge  of  50%  RH(0.628  mol/m3  at 
298  °K)  air  to  clean  bed  correspond  to  a  region  of  the  isotherm  that  is  unfavorable 
for  adsorption.  The  effluent  concentration  and  temperature  are  plotted  as  a 
function  of  time  in  Figure  4.  The  effluent  concentration  versus  time  curve  is  typical 
of  a  gradual  wave.  The  breakthrough  curve  is  concave  downward  to  the  time  axis. 
After  a  short  period  of  time,  an  appreciable  percentage  of  the  feed  concentration  is 
seen  in  the  effluent.  The  temperature  response  is  also  typical  of  a  gradual  wave. 

The  temperature  quickly  rises  to  a  maximum  (in  this  case  about  305  °K)  and  then 
falls  gradually  back  to  the  feed  temperature. 

Construction  of  solutions  using  equilibrium  theory  requires  the  use  of  the 
two  rules  stated  earlier.  A  wave  will  be  gradual  if  possible.  Slopes  in  the  q  versus 
v  c  plane  must  increase  from  initial  to  feed  conditions.  There  are  two  sets  of 
characteristics  corresponding  to  two  roots  of  the  equation  29.  The  two 
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characteristic  curves  are  required  in  order  to  obtain  the  path  that  connect  any 
initial  condition  to  any  feed  condition.  The  1  wave  is  the  first  wave  to  leave  the 
bed,  therefore  it  must  be  the  wave  connected  to  the  initial  condition  while  the  2 
wave  is  connected  to  the  feed. 

The  slope  of  the  1  characteristic  on  the  q  versus  v*c  plane  is  decreasing 
from  the  initial  condition,  a  clean  bed.  This  means  there  will  be  a  shock  path  that 
leaves  this  point.  This  is  shown  as  the  bottom  dotted  line  in  Figure  2  and  the 
upper  dotted  line  in  Figure  3  emanating  from  the  origin.  The  2  wave  that  connects 
the  feed  condition,  50%  RH,  has  a  constantly  increasing  slope  from  initial  to  feed 
condition  corresponding  to  a  gradual  wave  represented  by  the  bold  solid  line  in 
Figures  2  and  3.  The  point  where  the  shock  from  the  initial  condition  and  gradual 
wave  meet  can  be  found  by  stepping  out  in  temperature  along  each  path  using 
equations  29  and  30  until  they  intersect.  This  corresponds  to  a  plateau 
concentration.  Figures  2  and  3  make  it  easy  to  visualize  how  the  solution  path  cuts 
across  isotherms  producing  a  temperature  wave. 

There  is  very  close  agreement  between  the  stage  model,  the  equilibrium 
theory  model,  and  the  data  as  seen  in  Figure  4.  The  equilibrium  theory  solution 
does  predicts  a  larger  temperature  rise  than  the  stage  model.  This  is  expected 
because  equilibrium  theory  assumes  adiabatic  behavior.  The  actual  recorded 
temperature  does  not  reach  as  high  as  the  prediction  because  during  the  first 
minutes  of  the  run  the  sample  tube  was  periodically  taken  off  line  so  that  the  water 
uptake  could  be  recorded.  Note  that  both  the  concentration  and  the  temperature 
wave  travel  with  the  same  velocity  as  seen  in  Figure  4.  This  is  predicted  by 
equilibrium  theory  using  equation  28. 

A  preliminary  study,  using  only  isothermal  equilibrium  theory,  did  not  fit 
the  data  indicating  that  the  effect  of  temperature  is  very  important  to  this  process. 
There  is  a  gradual  wave  associated  with  the  low  RH  for  the  isothermal  case. 
However,  the  shock  associated  with  temperature  tends  to  overrun  these  gradual 
wave  characteristics  ill  the  adiabatic  case  leading  to  a  shock  in  temperature  and 
concentration.  Note  that  because  the  water  isotherm  is  unfavorable  at  lower  RH 
these  low  concentrations  tend  to  run  out  of  the  column  quickly.  The  large 
temperature  rise  associated  with  the  1  wave  has  the  effect  of  removing  water  from 
the  bed  that  was  just  adsorbed  at  the  lower  temperature.  This  increases  the 
required  time  for  complete  saturation  of  the  bed.  It  appears  that  isothermal 
adsorption  would  provide  a  faster  means  of  achieving  saturation. 
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t  Breakthrough  of  80%  RH  Feed  To  Bed  Equilibrated  At  50%  RH. 

The  water  isotherm  is  favorable  between  50%  RH(0.628  mol/m3  at  298°K) 
and  80%  RH(1.001  mol/m3  at  298°K)  as  seen  in  Figure  2.  This  is  the  typical  shape 
of  an  adsorption  isotherm,  and  one  could  expect  the  standard  S-shaped 
breakthrough  curve.  Figure  5  shows  the  results  of  an  80%  RH  feed  to  a  bed 
initially  equilibrated  at  50%  RH.  Although  the  shape  of  this  curve  is  not  very 
sharp,  the  breakthrough  behavior  does  follow  an  S-shaped  pattern. 

This  S-shaped  pattern  is  called  an  abrupt  wave  or  a  constant  pattern 
wave.  This  means  that  at  a  certain  distance  into  the  bed,  the  shape  of  the 
concentration  front  does  not  change  as  the  wave  is  transmitted  down  ihe  bed. 

There  is  also  an  initial  shock  characterized  primarily  by  a  temperature  rise. 

Referring  back  to  Figure  2,  the  2  characteristic  emanating  from  the  feed  of 
80%  RH  and  the  1  characteristic  from  the  initial  concentration  of  50%  RH  both 
violate  the  gradual  wave  criterion.  The  slopes  of  these  waves  decrease  from  initial  to 
feed  conditions,  therefore  each  of  these  waves  must  start  as  a  shock.  The 
equilibrium  theory  solution  then  requires  that  equation  30  be  solved  from  both  the 
initial  and  feed  conditions  until  they  intersect,  while  at  each  increment  checking  to 
see  whether  the  slope  of  the  gradual  wave  is  continuing  to  violate  the  gradual  wave 
criterion.  The  intersection  point  leads  to  a  plateau  condition.  This  path  is  close  to 
having  a  gradual  portion  because  it  is  so  close  to  the  inflection  point  of  the 
isotherm.  However,  the  adiabatic  nature  of  this  system  results  in  the  formation  of 
two  shocks. 

The  simulation  provided  by  the  stage  model  and  equilibrium  theory 
represent  the  data  well  for  this  run.  The  resulting  breakthrough  behavior  is  shown 
in  Figure  5.  Equilibrium  theory  provides  a  very  accurate  measure  of  the  center  of 
both  the  concentration  and  temperature  waves.  The  plateau  region  between  the 
two  shock  fronts  is  apparent  in  the  data.  Again  the  temperature  rise  measured  in 
the  experiment  is  reduced  due  to  sampling  of  the  uptake  weight  of  the  bed.  The 
effect  of  the  mass  and  heat  transfer  resistances  in  the  stageraodel  becomes  obvious 
in  the  shock  portion  of  the  breakthrough  curves.  These  resistances  tend  to  spread 
out  the  shock  wave  in  an  actual  experiment. 

4.4  Breakthrough  of  80%  RH  Feed  To  A  Clean  Bed. 

The  water  isotherm  changes  from  unfavorable  to  favorable  between  50% 
RH  and  80%  RH.  The  effect  of  this  inflection  in  the  isotherm  should  be  apparent 
on  the  breakthrough  behavior.  Figure  6  shows  the  breakthrough  data  measured  for 
80%  RH  feed  to  a  clean  bed.  Notice  that  the  resulting  breakthrough  curves  of  both 
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temperature  and  concentration  are  essentially  a  combination  of  the  first  two 
experiments.  The  only  difference  is  the  faster  flow  rate  (10  1/min  as  opposed  to  5 
1/min). 


The  equilibrium  theory  solution  path  is  characterized  by  an  initial  1  wave 
shock  as  in  the  run  with  50%  RH  feed  to  a  clean  bed,  the  lower  dotted  line  in 
Figure  3.  There  is  also  a  shock  emanating  from  the  feed  as  in  the  run  with  80%  RH 
feed  to  a  bed  equilibrated  to  50%  RH.  However  there  is  a  third  t  ransition.  As  the 
isotherm  changes  from  favorable  to  unfavorable  a  combined  wave  is  formed  that 
shows  both  shock  and  gradual  wave  behavior,  Figure  2.  It  is  possible  to  identify 
this  as  a  combined  transition  because  it  does  not  have  a  plateau  region  where  the 
two  waves  unite.  There  is  a  portion  of  the  gradual  wave  for  80%  RH  feed  to  clean 
bed  that  parallels  the  shock  path  for  80%  RH  feed  to  a  bed  equilibrated  at  50%RH. 
The  two  runs  are  not  identical  over  the  same  range  of  RH.  This  result  only 
becomes  apparent  in  the  adiabatic  analysis. 

The  shape  of  both  the  concentration  and  temperature  profiles  agrees  with 
the  stage  model  and  equilibrium  theory  predictions.  The  subtle  shape  of  the 
combined  transition  is  not  obvious  from  the  data.  Equilibrium  theory  allows  this 
combined  transition  to  be  identified. 

5.  CONCLUSIONS 

•  The  three  types  of  waves  that  can  be  formed  in  and  passed  through  an 
adsorption  bed  are  demonstrated  both  mathematically  and  experimentally. 

•  The  stage  model  seems  to  fairly  represent  the  breakthrough  using 
correlated  equilibrium  data  and  a  priori  calculated  rate  parameters. 

•  Equilibrium  theory  allows  the  exact  nature  of  the  transitions  to  be 
explored. 

•  The  inability  of  equilibrium  theory  to  be  used  with  heat  and  mass  transfer 
resistances  leads  to  some  discrepancy  between  this  method  and  the 
experimental  results. 
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NOTATION 


a 

A,  B,  C 

c* 

c 


Uf 

h 

s 

kv 

Kl,  K2 

L 

P 

q 


%at 

ri 

R 

t 

T 


To 


v 


V 

o 

z 

e 

C 

9 

A 


*b 

Pf 

p{o 


—  surface  to  volume  ratio(m"*) 

=  Antoine  equation  parameters  (p  in  mm  Hg) 

=  gas  phase  concentration  of  water(mol/m° ) 

=  gas  phase  concentration  of  water  in  the  pellet(mol/m  ) 
=  heat  capacity  of  adsorbed  phase (kJ/kgK) 

=  heat  capacity  of  gas  phase(kJ/kgK) 

=  heat  capacity  of  solid  phase  (kJ/kgK) 

=  enthalpy  of  gas  phase(kJ/kg) 

=  enthalpy  of  solid  phase(kJ /kg) 

=  mass  transfer  coefficient  (m/s  Ref.  16) 

=  coefficients  of  water  adsorption  isotherm,  Eq.  (32) 

=  length  of  bed(m) 

=  partial  pressure  of  water(Torr) 

=  solid  phase  concentration  (mol/kg) 

=  solid  phase  concentration  at  saturation(mol/kg) 

=  inner  radius  of  filter  tube(m) 

=  gas  constant 
=  time 

=  temperature(K) 

=  temperature  of  feed  gas(K) 

=  heat  capacity  reference  temperature  (K) 

=  overall  heat  transfer  coefficient  of  filter  tube(kJ/m  Ks) 
=  interstitial  velocity(m/s) 

=  interstitial  velocity  of  feed(m/s) 

=  axial  distance  coordinate  (m) 

=  void  fraction  of  packing 
=  dimensionless  axial  coordinate,  Eq.  (7) 

=  fractional  saturation  of  the  bed 
=  heat  of  adsorption(kJ /mol) 

=  bulk  density  of  packing(kg^m^) 

=  density  of  gas  phase  (kg/m0) 

=  density  of  inlet  gas  (kg/m^) 

=  dimensionless  time,  Eq.  (6) 
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Figure  1.  Schematic  of  Breakthrough  Apparatus 
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Figure  2.  Path  In  Modified  Isotherm  Plane 

50ZRH,  80 '/.RH  Feed  To  Clean  Bed 
And  80*RH  Feed  To  5C£RH  Bed  Of  BPL 


.  0008  r 


Blank 


28 


LITERATURE  CITED 


1.  Doong,S.J.,  and  Yang,R.T.,  'Adsorption  of  Mixtures  of  Water  Vapor  and 
Hydrocarbons  By  Activated  Carbon  Beds:  Thermodynamic  Model  for  Adsorption 
Equilibrium  and  Adsorber  Dynamics,'  AIChE  Symposium  Series:  Adsorption  and 
Ion  Exchange  vol.  259(83),  p.  87  (1987). 

2.  Rodrigues,  A.E.  and  Beira,  E.C.,  'Staged  Approach  of  Percolation 
Processes,'  AIChE  J.  vol.  25,  p.  416  (1976). 

3.  Ikeda,  X.,  'Performance  of  the  Nonisothermal  Fixed-Bed  Adsorption 
Column  with  Nonlinear  Isotherms,'  Chem.  Eng.  Set.  vol.  34,  p.  941  (1981). 

4.  Friday,  D.K.,  and  LeVan,  D.,  'Solute  Condensation  in  Adsorption  Beds 
During  Thermal  Regeneration,'  AIChE.  J.  vol.  28,  p.  86  (1982). 

5.  Friday,  D.K.,  and  LeVan,  D.,  'Thermal  Regeneration  cf  Adsorption  Beds: 
Equilibrium  Theory  for  Solute  Condensation,'  AIChE  J.  vol.  30,  p.  679  (1984). 

6.  Acrivos,  A.,  'Method  of  Characteristics  Technique:  Application  to  Heat 
and  Mass  Transfer  Problems,'  Ind.  Eng.  Chem.  vol.  48,  p.  703  (1956). 

7.  Pan,  C.  Y.,  and  Basmadjian,  D.,  "Constant  Pattern  Adiabatic  Fixed-Bed 
Adsorption,"  Chem.  Eng.  Set.  vol.  22,  p.  285  (1967). 

8.  Pan,  C.Y.,  and  Basmadjian,  D.,  "An  Analysis  of  Adiabatic  Sorption  of 
Single  Solutes  in  Fixed  Beds:  Equilibrium  Theory,'  Chem.  Eng.  Sci.  vol.  26,  p.  45 
(1971). 

9.  Rhee,  H.K.,  Aris,  R.  and  Amundson,  N.R.,  First  Order  Partial  Differential 
Equations:  Volume  I,  Theory  and  Applications  of  Single  Equations,  pp.  188-410, 
Prentice-Hall,  Englewood  Cliffs,  NJ,  1986. 

10.  Ruthven,  D.M.,  Principles  of  Adsorption  and  Adsorption  Processes, 
pp.  274-323,  John  Wiley  and  Sons,  New  York,  NY  1984. 

11.  De Vault,  D.,  "The  Theory  of  Chromatography,'  J.  Amer.  Chem.  Soc. 
vol.  65,  p.  532  (1943). 

12.  Lax,  P.D.,  'The  Formation  and  Decay  of  Shock  Waves,'  Amer.  Math. 
Monthly  vol.  79,  p.  227  (1972). 


29 


13.  Hacskaylo,  J.J.,  and  LeVan,  M.D.,  'Correlation  of  Adsorption  Equilibrium 
Data  Using  a  Modified  Antoine  Equation:  A  New  Approach  for  Pore-Filling 
Models,'  Langmuir  vol.  1,  p.  97  (1987). 

14.  Reid,  R.C.,  Prausnitz,  J.M.,  and  Sherwood, T.K.,  The  Properties  of  Gases 
And  Liquids ,  3rd  ed.,  p.  634,  McGraw-Hill,  New  York,  NY,  1977. 

15.  Mahle,  J.J.,  and  Friday,  D.K.,  Water  Adsorption  Equilibria  On 
Mieroporous  Carbons ,  CRDEC-TR-018,  U.S.  Army  Chemical  Research, 
Development  and  Engineering  Center,  Aberdeen  Proving  Ground,  MD,  November 
1988,  UNCLASSIFIED  REPORT. 

16.  Ranz,  W.E.  and  Marshall,  W.R.,  'Evaporation  From  Drops  Part  II,' 
Chem.  Eng.  Prog.  vol.  48,  p.  173  (1952). 


